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Abstract To clarify the structure and organization of butterfly communities, I analyzed a 
transect-count data set of seven butterfly communities in varying types of woodland habitat at 
the northern foot of Mt Fuji, central Japan. The results of statistical analyses including 
UPGMA suggested that the butterfly community was composed of four species subgroups that 
showed either relatively high abundances and wide local distributions (subgroup IV), relatively 
low abundances and wide local distributions (subgroups I and IID, or relatively low abundances 
and restricted local distributions (subgroup II). The result of PCA suggested that the extent of 
growth in the mantle and sleeve plant communties (i e. the degree of development of forest edges 
(corresponding to ecotone)) in each habitat played the major role in the diversification, forma- 
tion, and organization of the four species subgroups. The two species subgroups IV and II, 
although quite different in abundance and distribution patterns, both had a great influence on 
the variations in the community indices such as species richness (total number of species) and 
total population abundance among the seven transect sites, whereas the other two subgroups I 
and III did not. From the analysis of the habitat types associated with the component species 
of the four species subgroups, it may be suggested that the differences in population abundances 
among the four species subgroups could be interpreted in terms of the sources and sinks theory 
(e. g. Pulliam, 1988). Overall, my analysis demonstrates that statistical approaches, including 
multivariate analysis, are very effective in the analysis of the structure and organization of 
butterfly communities. 
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Introduction 


Up to the present, there have been many studies on the diversity, structure, and composition 
of butterfly communities (see Gilbert & Singer, 1975; Gilbert, 1984; Yamamoto, 1988; Porter 
et al., 1992; Ishii, 1993; Pollard & Yates, 1993; New, 1997; Imai & Ishii, 1998). In these 
studies, butterfly communities have often been analyzed from the viewpoint of community 
parameters such as species richness, species diversity, equitability, total population abundance, 
etc. (see New, 1997; Natuhara, 1998; Sakuratani, 1998) and ‘environmental indices’ (e. g. 
Tanaka, 1988; Sunose, 1993, 1996, 1998; Hattori et al., 1997; Tashita & Ichimura, 1997). 
However, especially in Japan, few attempts have been made to analyze butterfly community 
structure using a multivariate statistical approach, except by Fujii (1988), Yamamoto (1991a, 
1991 b, 1993, 1996) and Natuhara (1996). In addition, there have been few studies that have 
tried to determine the organizing force of butterfly communities using statistical approaches. 


In my previous study (Kitahara, 1999), I monitored seven butterfly communities in varying 
types of woodland habitat by the line transect method during 1996 at the northern foot of Mt 
Fuji, central Japan. In this study, I used the butterfly community data set stated above, and 
analyzed the structure and organization of the communities using a multivariate statistical 
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approach. My goals of this study were 1) to discriminate the species subgroups which 
constituted the butterfly community, 2) to identify the environmental factors contributing to 
the diversification and formation of the species subgroups of the community (if they are 
discovered), and 3) to clarify the characteristics and role in a community of the species 
subgroups (if they are discovered). 


Materials and methods 


Data set analyzed 


In this paper, I used the data collected in my previous study; because details of the census sites 
are described elsewhere (Kitahara, in press), only short descriptions follow. I established 
seven butterfly transects in different woodland types (alt. ca 1,000 m) at the northern foot of 
Mt Fuji (alt. 3,776 m), Yamanashi Prefecture, central Japan. The two sites at Tsuchimarubi 
were located along the outer (site 1) and inner (site 2) forest edges of broad-leaved deciduous 
plantations that were sunny and relatively undisturbed. The Nakanochaya site was the most 
shady transect, located in a coniferous and deciduous mixed forest. Тһе Hokuroku-Park site 
was in an area of naturally remaining and artificially transplanted woods with human- 
managed undergrowth. The Kankyouken site was in a recently logged and highly disturbed 
area surrounded by a secondary red pine forest. The two sites at Funatsutainai were located 
along the outer (site 1) and inner (site 2) forest edges of coniferous plantations that were sunny 
and relatively undisturbed. Table 1 summarizes the characteristics of the seven transect sites. 


I monitored butterfly communities along each transect twice a month during the adult flight 
season (from April to November, 1996) by transect counts (Pollard, 1977; Yamamoto, 1988; 
Ishii, 1993; Pollard and Yates, 1993). Details of the census method are described elsewhere 
(Kitahara, 1999). I observed a total of 72 butterfly species along the seven transects. Pieris 
melete and Pieris napi were treated as a single species (Pieris sp.) because of the difficulty in 
distinguishing them in field censuses, especially during flight. 


Data analysis 


I analyzed butterfly community structure based on the population abundance of each species 
along each transect. Population abundance of each species was expressed as the mean 
monthly density for each transect, which was obtained by averaging the monthly densities, 
only when butterfly species were found. The monthly density was usually the mean of the 
two observations in each month. Exclusion of months when no butterflies were observed 
minimized the effect of variations in the voltinism of butterfly species on the yearly abundance 
estimate (i е. the mean monthly density). Average population abundance of each compo- 
nent species in the community was calculated by averaging the population abundances of 
each transect, only where butterfly species were recorded, and converting this value into that 
per 1 km. Exclusion of transects where no butterflies were recorded minimized the effect of 
different distribution patterns (widespread or restricted) of butterfly species on the average 
population abundance. All butterfly species observed, the population abundance of each 
species at each transect, and the average population abundance of each species in the entire 
butterfly community are described in the Appendix. 


A 72 (species) by 7 (transect sites) matrix based on the population abundance (mean monthly 
density) of each species at each transect was subjected to multivariate analyses. Species 
composition and clustering were examined by cluster analysis (unweighted pair-group 
method using arithmetic average (UPGMA)) based on a 72 by 72 Euclidean distances matrix. 
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Study site 





Altitude 
Vegetation 
High tree layer 


Subhigh tree layer 


Tsuchimarubi | 
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Summary description of the seven study sites. 





Tsuchimarubi 2 





Nakanochyaya 





920-940 m 


Alnus hirsuta plantation 














920-940 m 


Alnus hirsuta plantation 


1,050-1,075 m 


Pinus densiflora & Larix 
leptolepis plantation 


Broad-leaved deciduous 
trees 























Degree of growth of Intermediate High Intermediate 
mantle community 
Degree of growth of High High Low 
sleeve community 
Sunlight condition (forest More bright Bright Highly dark 
edge) 
Level of human Intermediate Low Low 
disturbance (forest edge) 
Hokuroku park Kankyoken Funatsutainai 1 Funatsutainai 2 
1,010-1,030 m 1,030-1,050 m 1,050- 1,060 m 1,050-1,070 m 


Remaining or transplanted 
trees 


Remaining or transplanted 
trees 

Low 

Low 


Bright 


Very high 


Pinus densiflora secondary 
forest 


Broad-leaved deciduous 
trees 

Very low 

Very low 

Bright 


Very high 


Pinus densiflora & Larix 
leptolepis plantation 


Broad-leaved deciduous 
trees 

High 

Very high 

Bright 


Intermediate 











Pinus densiflora & Larix 
leptolepis plantation 


Broad-leaved deciduous 
trees 

High 

High 

Bright 


Low 





Species ordination was carried out by principal components analysis (PCA) based on a 
variance-covariance matrix. Characteristics of the species groups or clusters found through 
multivariate approaches were examined in relation to biological attributes of the component 
species. Those used in the analyses were 1) voltinism (based on Unno and Aoyama (1981), 
and also on the actual data of seasonal changes in the observed number of individuals of each 
species in my study), 2) larval diet breadth (based on Endo and Nihira (1990)), 3) overwinter- 
ing stage (based on Unno and Aoyama (1981)), and 4) habitat type (based on Miyatake 
(1998)). 


Results 


Species composition and grouping by cluster analysis (UPGMA) 


Figure 1 shows the community dendrogram of Euclidean distances between all butterfly 
species. When the level of the Euclidean distance for the division of samples was determined 
arbitrarily to be 2.55, there were four major clusters (species subgroups) (I, II, Ш, and other 
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species (IV)). The mean values of average population abundance of component species 
among the four species subgroups were significantly different (ANOVA: F—29271, P<0.01) 
(Fig. 2a). In addition, the mean values for the numbers of transect sites where each compo- 
nent species was observed were also significantly different among the four species subgroups 
(ANOVA: F=23.768, Р< 0.01) (Fig. 2b). 


These results implies that the butterfly community was composed of the following 4 species 
subgroups with different densities and/or distributions: 1) a species subgroup with high 
abundances and wide distributions (IV), 2) that with low abundances and wide distributions 
(I and Ш), and 3) that with low abundances and narrow distributions (II). 


Species ordination by PCA 


Distribution of the butterfly species observed on major (the first)- and minor (the second)-axes 
plane by PCA based on a variance-covariance matrix is shown in Fig.3. The species 
subgroups (I, П, Ш and IV) detected by cluster analysis distributed from left to right sides 
along the major axis. The tightly packed species subgroup II distributed mainly on the left 
side of the major axis. On the contrary, the scattered species subgroup IV distributed mainly 
on the right side of the axis. The other two subgroups I and III were almost situated between 
subgroups II and IV along the major axis. These PCA results support the validity of the 
species grouping and clustering detected by UPGMA stated above. 


The cumulative contribution of the first and second principal components was 74.05% (Table 
2). For the first axis, the values of eigenvectors were all positive. The Funatsutainai sites 
1 and 2, where more growth mantle and sleeve plant communities were featured (Table 1), 
had the highest eigenvector values (>0.5), while the Kankyouken site, which was an area 
with almost no mantle and sleeve plant communities (because it was recently logged and 
highly disturbed by humans) (Table 1), had the lowest eigenvector value (0.07). The degree 
of growth in the mantle and sleeve communities in the other four transect sites was almost 
intermediate between the Funatsutainai two sites and the Kankyouken site. These results 
suggest that the first axis reflected a gradient of the extent of growth in the mantle and sleeve 
plant communities of the sites. 


Table 2. Eigenvalues and eigenvectors of the variance-covariance matrix used for PCA of the 
butterfly community. 








Principal components 


Axis 1 Axis 2 














Eigenvalue 52.2230 12.7664 
Contribution to community variance (%) 59.50 14.55 
Cumulative 96 59.50 74.05 
Eigenvector Tsuchimarubi 1 0.2727 0.1443 
Tsuchimarubi 2 0.3917 0.7331 
Nakanochaya 0.2778 —0.5527 
Hokuroku-Park 0.2259 — 0.0389 
Kankyoken 0.0725 — 0.0730 
Funatsutainai 1 0.5882 0.0567 
Funatsutainai 2 0.5411 — 0.3552 
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Fig. 1. Dendrogram of 72 butterfly species classified by the cluster analysis (unweighted 
pair-group method using arithmetic average (UPGMA)) based on a 72 by 72 Euclidean 
distance matrix. The species numbers in the figure correspond to those in Appendix. 


For the second axis, the first and second highest eigenvector values were detected in 
Tsuchimarubi sites 2 and 1, which were located in the inner and outer areas of a broad-leaved 
deciduous forest, respectively (Table 1). On the other hand, the first and second lowest 
values were found in the Nakanochaya site and the Funatsutainai site 2, which were located 
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Fig. 2. 


Average population abundance 


Number of sites observed 


Masahiko KITAHARA 






50 
B Mean 


О Mean+o 
X Меап– о 








40 


30 


20 


10 


IV lll | [ 
Species subgroup 


ШП Mean 
O Mean+o 
X Mean- 






IV it | | 
Species subgroup 


Differences in the average population abundance (a) and the number of census sites 
observed (b) among the four species subgroups of the butterfly community. 


in the inner and outer areas of a conifer-dominated forest, respectively (Table 1). These 
results suggest that the second axis reflects a gradient from broad-leaved deciduous woods via 
the ecotonal areas of woods to conifer-dominated woods. 


The results of PCA indicate that the species grouping and clustering by UPGMA are 
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Fig. 3. Scattergram of 72 butterfly species along the first and the second principal component 
axis based on PCA ordination. The species numbers in the figure correspond to those 
in Fig. 1 and Appendix. 


governed by the first axis (the extent of growth in the mantle and sleeve plant communities 
of the sites). 


Correlations between butterfly species richness and abundance in each transect site, and 
the number of species in each subgroup in each transect site 


Butterfly species richness (the total number of species in each transect site) was significantly 
and positively correlated with the number of species in subgroups II and IV in each transect 
site (П: 7—0.988, P<0.01; IV: r=0.881, Р<0.01), but not with that in subgroups I and III 
(I: 7—0.003, P>0.05; Ш: r=0.735, P>0.05) (Fig. 4а). In addition, butterfly species abun- 
dance (the total population abundance of all species observed in each transect site) was 
significantly and positively correlated with the number of species in subgroups П and IV in 
each transect site (II: 7—0.894, P<0.01; IV: r=0.963, Р< 0.01), but not with that in sub- 
groups I and Ш (E: r= —0.372, P>0.05; Ш: r=0.615, P>0.05) (Fig. 4b). 


Comparisons of various characteristics among the four butterfly subgroups 


Table 3 shows the comparisons of various ecological characteristics cited from the literature 
among the 4 butterfly subgroups detected in the previous analyses. As for the voltinism of 
the component species, subgroup III was composed completely of multivoltine species. In 
contrast, subgroup IV was composed mostly (81.896) of univoltine species. The other two 
subgroups I and II were composed of multi- and univoltine species in similar proportions 
(about 50.0% each). There were no marked differences in the overwintering stage of 
component species among the four subgroups; the highest values for the overwintering stage 
was obtained in the larval stage for all subgroups I (43.8%), П (37.5%), III (60.0%), and IV 
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Fig. 4. Relationships of the number of species in each subgroup with species richness (the total 
number of species) (a), and with total population abundance (b). 


(59.1%). 


Among the categories of habitat type as appearing in Table 3, the habitat type usually 
occurring in the present study area is considered to be broad-leaved deciduous warm- 
temperate forest, broad-leaved deciduous cold-temperate forest, dry-type grassland, and sparse 
woods (marked by ** in Table 3. When we note the total proportion of the habitat types 
occurring in the present study area as stated above (see Table 3), 81.6% of well-known habitat 
types associated with subgroup IV as cited from the literature corresponded to the habitat 
types occurring in the present study area, whereas the well-known habitat types associated 
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Table 3. Comparisons of various characteristics cited from the literature among the four 
species subgroups of the butterfiy community. 


Species subgroup 





I II Ш IV 

Characteristics No. spp. (%) No. spp. (%) No. spp. (%) No. spp. (%) 
Voltinism & larval diet breadth* 

Multivoltine & generalist 2 (25.0) 9 (18.8) 1 (20.0) 1 (9.1) 

Multivoltine & specialist 2 (25.0) 14 (29.2) 4 (80.0) 1 (9.1) 

Univoltine & generalist 0 (0.0) 6 (12.5) 0 (0.0) 4 (36.4) 

Univoltine & specialist 4 (50.0) 19 (39.6) 0 (0.0) 5 (45.5) 
Overwintering stage 

Egg 0.5 (6.3) 11 (22.9) 0 (0.0) 25 (227) 

Larva 35 (43.8) 18 (37.5) 3 (60.0 | 6.5 (591) 

Pupa 3 (375) n (229) 1  Q00 1 (9.1) 

Adult 1 (12.5) 7 (14.6) 1 (20.0) 1 (9.1) 

Unknown 0 (0.0) 1 (2.1) 0 (0.0) 0 (0.0) 
Habitat type 

Broad-leaved evergreen forest 2.5 (31.3) 7 (14.6) 1 (20.0) 0.5 (4.6) 

Broad-leaved deciduous warm-temperate 2 (25.0) 73 (15.3) 0 (0.0) 3 (27.3) 

forest** 
Broad-leaved deciduous cold-temperate 0 (0.0) 43 (9.0) 0 (0.0) 1.5 (13.6) 
forest** 
Upper river basin 1 (12.5) 6 (12.5) 0 (0.0) 0 (0.0) 
Lower or middle river basin 0.5 (6.3) 3.3 (6.9) 2 (400) 0.5 (4.6) 


Dry-type grassland** 2 (25.0) 5.5 (11.5) 1.5 (30.0) 4 (36.4) 
Wet-type grassland 0 (0.0) 1.5 (3.1) 0 (0.0) 0 (0.0) 
Sparse wood** 0 (0.0) 3.5 (7.3) 0 (0.0) 0.5 (4.6) 
Arable field 0 (0.0) 6 (12.5) 05 (100) 0.5 (4.6) 
0 
0 


Residential and urban area (0.0) 1 (2.1) 0 (0.0) 0.5 (4.6) 


Total proportion of habitat types (50.0%) (43.0%) (30.0%) (81.6%) 
occurring in the study area 
(marked by**) 


Total proportion of habitat types (50.0%) (57.0%) (70.0%) (18.4%) 
rarely occurring in the study area 


* Larval diet breadth: specialist; the species whose larvae had been reported to feed on hostplants 
belonging to one taxonomic family, generalist; the species whose larvae had been reported to feed on 
hostplants belonging to two or more taxonomic families. 

**: the habitat type occurring in the study area. 


with subgroup III such as lower and middle river basin and arable fields rarely occurred in 
the present study area (only 30.0% occurred). For the other subgroups I and II, about half 
of the known habitat types cited from the literature occurred regularly in the study area, while 
the other half such as broad-leaved evergreen forest and upper river basin did not. 
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Discussion 


Species composition of the butterfly community 


Results of the cluster analysis (UPGMA) and ANOVAs of average population abundances 
and local distribution patterns of component species suggested that the butterfly community 
was composed of four species subgroups. These subgroups were classified into three types: 
1) that with high abundances and wide distributions (one subgroup), 2) that with low 
abundances and wide distributions (two subgroups), and 3) that with low abundances and 
narrow distributions (one subgroup). In addition, the results of PCA also roughly support- 
ed the species grouping and clustering generated by UPGMA. The 4 species subgroups 
discriminated by UPGMA were arranged almost independent of each other along the first 
axis of the PCA. This suggests that the extent of growth in the mantle and sleeve plant 
communties (i e. the degree of development of forest edges (corresponding to ecotones)), 
which related to the first axis of the PCA, played the central role in the diversification, 
formation, and organization of the species subgroups of the butterfly community. 


Up to the present, it has been generally thought that the development of ecotones (forest edges 
in woodlands) enhances species diversity and richness in communities (e. g. Pianka, 1988; 
Smith, 1990; Ishii et al., 1995; Ishii, 1996; Maeda, 1996; Natuhara, 1996; Tanaka & Arita, 
1996; Yata, 1996). This is because such ecotonal communities include various characteristic 
types of species as a result of a variety of adjacent habitat components (Pianka, 1988; Smith, 
1990; Yamamoto, 1993; Maeda, 1996). For butterflies, it has been pointed out that 
increases in resources such as nectar plants for adults and hostplants for larvae, based on the 
development of ecotones, enhances species diversity and richness in ecotonal butterfly 
communities (Ishii et al., 1995, Ishii, 1996; Kitahara, unpublished). In addition to these 
mechanisms, my analysis demonstrates that the growth of mantle and sleeve plant com- 
munties, based on the development of forest edges, promotes the diversification and formation 
of species groups in butterfly communities, and as a result, enhances the diversity and species 
richness of ecotonal butterfly communities. This mechanism has received little attention so 
far. In general, it is considered that the growth of mantle and sleeve plant communities is 
governed mainly by physical disturbance and stress (e. g. Grime, 1977; Southwood, 1988). 
Therefore, it is thought that the physical conditions of disturbance and stress also govern the 
diversity and species richness of butterfly communities. Indeed, many studies have shown 
that butterfly species diversity is largely influenced by the degree of disturbance (e. g. Spitzer 
et al., 1993, 1997; Kitahara & Fujii, 1994; Blair & Launer, 1997), and of stress (sunlight 
condition) (е. g. Warren, 1985; Greatorex-Davies et al., 1993; Sparks et al., 1996). 


On the other hand, my analysis showed that the differences in the component tree species of 
woodlands, reflecting the second axis of PCA, played little role in the diversification and 
formation of the species subgroups of the butterfly community. In birds, it is well known 
that species diversity is more affected by foliage height diversity (indicating the extent of the 
development of stratification) than plant species diversity (е. g. MacArthur & MacArthur, 
1961; MacArthur et al., 1966; Cody, 1970; Karr, 1971; Higuchi,1978). Butterflies, on the 
other hand, in both the adult and larval stages depend mostly on plants as their necessary 
resource, and therefore, butterfly species diversity is generally thought to be strongly 
influenced by plant species diversity (as suggested by Katano (1995) and Ishii (1996)). My 
analysis, however, showed that the diversification and formation of species subgroups in the 
butterfly community, which contributed to the enhancement of species diversity, was promot- 
ed by the development of forest edge plant communties rather than the distinction of forest 
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component tree species. The fact that this phenomenon was found in herbivorous insects 
such as butterflies is important. Probably, this is due to the high dependence of the nectar 
plants and hostplants of butterflies on the sleeve and mantle communities in forest edges. 
Similar arguments have been made by Ishii et al. (1995). 


Consequently, the butterfly community was found to be composed of several species sub- 
groups with different characteristics, and it was concluded that the degree of development in 
the habitat structure of forest edges may play an important role in the diversification and 
formation of species subgroups. 


Characteristics of the species subgroups of the butterfly community 


Table 4 summarizes the characteristics of the four species subgroups of the butterfly commu- 
nity based on the analyses of this study. As discussed above, it is suggested that these 
subgroups diverged along a gradient of the degree of development of forest edge sleeve and 
mantle plant communties. The species subgroup IV at one end of the gradient was 
associated with highly developed forest edge plant communities, and had large population 
sizes and wide local distributions. In contrast, subgroup II at the other end of the gradient 
was probably associated with components other than forest edge communities such as forest, 
grass, arable fields, and urban areas, and had small population sizes and restricted local 
distributions. The other two subgroups І and Ш were positioned between the two sub- 
groups stated above along the gradient, and had intermediate features (small population sizes 
and wide distributions). 


Table 4. Summary of the characteristics of the four species subgroups found in the butterfly 
community at the northen foot of Mt Fuji, central Japan. 











Species subgroup 

















Characteristics I II Ш IV 
Population density Low Low Low High 
Local distribution pattern Wide Narrow Wide Wide 
PCA I (The degree of develop- Intermediate Low Intermediate High 
ment of mantle and sleeve 
communities)” 
PCA II (Broad-leaved wood- co- Broad-leaved? Broad-leaved? Broad-leaved? Broad-leaved? 
niferous wood gradient)? Coniferous Coniferous Coniferous Coniferous 


Correlation of the no. of spp. in 
each group? 


with species richness NS seat NS ane 
with total population abun- NS ++ NS ++ 
dance 
Voltinism? — — M (U) 
Occurrence rate of the main hab- 50.0% 43.0% 30.0% 81.6% 


itats of each species subgroup in 
the study area 











1) First principal component. 

2) Second principal component. 

3) Subgroup included both broad-leaved and coniferous woody areas depended species. 

4) NS: not significant, ++: P<0.01 

5) M: АП species were multivoltine, (U): Most species were univoltine. —: No trends were detected. 
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According to the results of the correlation analysis between butterfly species richness and 
abundance in each transect site, and the number of species in each subgroup in each transect 
site, only the number of species in subgroups II and IV were significantly and positively 
correlated with butterfly species richness (total number of species) and abundance (total 
population abundance of all species) in each transect site. This reflects the increase in both 
butterfly species richness and abundance with an increase in the number of species in 
subgroups П and IV. The presence of two types of species groups, those contributing to and 
those irrelevant to community indices such as species richness and total population abun- 
dance in a butterfly community, has been reported elsewhere (e. g. Kitahara & Fujii, 1994). 
Moreover, it is noteworthy that the subgroups II and IV, with quite distinct features from each 
other, played the major roles in the determination of certain kinds of community indices of 
the butterfly community. Probably, this is due to the results which show that the differences 
in habitat structure among the study sites highly reflected both the differences in forest edge 
habitat structure and those in non-forest edge habitat structure such as woodlands, grasslands, 
and plain fields among the study sites. The Tsuchimarubi sites, which showed the greatest 
species-richness and population-abundance in this study, аге thought to have provided 
suitable habitat structure for many butterfly species of the subgroups IV , associated with 
developed forest edges, and II, more associated with non-forest edge components such as 
forest, grass, and plain fields. This is consistent with the arguments that greater species 
richness in forest edges is a result of the presence of various types of species groups characteris- 
tic to forest edge and characteristic to the neighboring habitats such as woodlands, grasslands, 
and plain fields (Pianka, 1988; Smith, 1990; Maeda, 1996). 


The main habitat types of butterflies usually occurring in the study area were more consistent 
with those of the subgroup IV based on the literature than those of the other 3 subgroups. 
This means that, for most species of subgroup IV, the present study area may have been 
located near the center of their geographic ranges (? e. the source area). On the other hand, 
there is a possibility that, for most of the species in the other 3 subgroups, the present area 
corresponded to the periphery of their ranges (i e. near the sink area). Usually, it is thought 
that source populations in source habitats produce an excess of individuals (population 
growth rate: r>0), whereas sink populations in sink habitats have less local productivity 
than local mortality (r« 0) (е. g. Pulliam, 1988; Lawton et al., 1994; Lawton, 1995; Pulliam 
& Dunning, 1997). In addition, the intrinsic rate of natural increase (r) is usually propor- 
tional to population abundance. Thus, the differences in the average population abundance 
among the 4 species subgroups were probably due to differences in the population growth rate 
(r) depending on whether the individuals were part of source or sink populations. That is, 
the subgroup IV with high abundances may have been associated with a higher r because it 
was part of a source population, whereas the other 3 subgroups with low abundances may 
have been associated with a lower r because they were part of sink populations. 


Consequently, in the butterfly community, there were different types of species subgroups, 
some of which played a great role in the determination of certain kinds of community indices, 
and the differences in the population abundances among the species subgroups may be 
interpreted in terms of the sources and sinks theory (e. g. Pulliam, 1988), as discussed above. 
Overall, my analysis demonstrates that statistical approaches, including multivariate analysis, 
is very effective in the analysis of the structure and organization of butterfly communities. 
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